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Abstract. Numerical calculations of lattice thermal conductivity are reported 
for the binary alloys NiPd and NiPt. The present work is a continuation of 
an earlier paper by us [1] which had developed a theoretical framework for the 
calculation of configuration-averaged lattice thermal conductivity and thermal 
diffusivity in disordered alloys. The formulation was based on the augmented 
space theorem [2] combined with a scattering diagram technique. In this paper we 
shall show dependence of the lattice thermal conductivity on a series of variables 
like phonon frequency, temperature and alloy composition. The temperature 
dependence of k(T) and its relation to the measured thermal conductivity is 
discussed. The concentration dependence of k appears to justify the notion of a 
minimum thermal conductivity as discussed by Kittel, Slack and others [3, 4], We 
also study the frequency and composition dependence of the thermal diffusivity 
averaged over modes. A numerical estimate of this quantity gives an idea about 
the location of mobility edge and the fraction of states in the frequency spectrum 
which is delocalized. 



PACS numbers: 72.15.Eb, 66.30.Xj, 63.50.+X 

1. Introduction 

Lattice thermal conductivity of substitutionally disordered alloys yield valuable 
information about the interactions of thermal excitations, phonons, with composition 
fluctuations on their crystal lattice. Over the past few years numerous experimental 
studies [5]- [8] of the thermal conductivity of disordered alloys have provided 
considerable insight into the nature of their elementary excitations. 

The theory of lattice thermal conductivity for perfect crystals and ordered alloys 
has been set up on a rigorous basis. However the same is not true for disordered alloys. 
The presence of disorder results in scattering that not only depends on the impurity 
concentration but also crucially on both the relative masses and size difference between 
the constituent atoms. For large mass or size differences, the effect of disorder can be 
quite unusual. Because of this, detailed comparison between theory and experiment 
on the basis of realistic models has not been very extensive. Model calculations are 
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mostly based on mean-field approaches with diagonal disorder alone. Whereas in 
phonon problems essential off-diagonal disorder in the force constants cannot be dealt 
within single site mean field approximations. Such disorder effects cannot be ignored 
in realistic calculations. 

When we come to comparison with experiment, however, we face a different kind 
of difficulty. In any experiment, the measured thermal conductivity k consists of the 
sum of an electronic component K e and a lattice component kl and k = n e + kl- 
Assuming the thermal analogue of Matthiessen's rule to be valid, the electronic 
thermal resistivity W e = l/n e is given by the sum of an ideal resistivity and an 
impurity or residual term. It is often assumed that the ideal resistivity remains 
unaltered by alloying and can be obtained from the measurements on pure metals. 
The residual resistivity W r can be calculated with the help of Widemann-Franz law 
W r — po(T)/(L T), where L is the Lorenz number and po(T) is determined by 
measuring the electrical resistivity of the alloy at several temperatures. The lattice 
component kl can then be separated out from the observed conductivity n. Overall 
what we would like to convey from these details is that a direct measurement of the 
lattice component of the thermal conductivity is not feasible. There always exists 
certain assumptions behind the calculation of n e and hence it is not possible always to 
obtain reliable estimates for this quantity which consequently affects the separation 
of Kl from the observed conductivity k. 

In an earlier communication [1] (hereafter referred to as AM) we had developed 
a formalism for the calculation of configuration averaged lattice thermal conductivity 
and thermal diffusivity for a disordered binary alloy. Unlike single-site mean-field 
approaches, this formulation has no restrictions on what kind of substitutional disorder 
can be studied. It explicitly takes into account fluctuations in masses, force constants 
and heat currents between different nuclei. It also maintains, on the average, the 
sum rule between the diagonal and off-diagonal parts of the dynamical matrix which 
eliminates the translational mode. We had shown that the dominant effect of disorder 
is to renormalize each propagator as well as the current terms in the Kubo formula. 

The purpose of this article is to implement the formulation (AM) for a detailed 
numerical study of the configuration averaged lattice conductivity and thermal 
diffusivity for disordered NiPt and NiPd alloys. The choice of the alloy systems is 
not arbitrary. We have chosen these alloys in order to propose a systematic study of 
the effects of mass disorder and strong force constant disorder. In the NiPd alloys, 
mass disorder is much larger than the force constant disorder. However both kinds of 
disorder dominate in the NiPt alloys. Moreover in both the two alloy systems there is a 
large size mismatch between the constituents. This indicates that the standard single 
site mean field theories would be inadequate to capture these effects. These alloy 
systems are therefore ideal for the illustration of the advantages of the augmented 
space block recursion method proposed by us [9]. Theoretically there have been 
many attempts [10]- [13] to develop an adequate approximation for understanding 
the lattice thermal conductivity of metals and alloys. The majority were based on 
model calculations either for perfect crystals or ordered alloys. One of the most 
successful mean-field approximation was the coherent potential approximation (CPA) 
[14] in combination with the appropriate Kubo formula. The CPA is a single site 
mean-field theory capable of dealing only with mass disorder. There have been many 
attempts to generalize the CPA in order to treat both off-diagonal and environmental 
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disorder. In many cases these approximations did not satisfy the necessary herglotz 
analytic and lattice translational invariance properties of the configuration averaged 
Green functions. Two very similar approaches, both based on the augmented space 
theorem, have been recently proposed : the itinerant cluster CPA (ICPA) [15] and 
the augmented space recursion [16]. These have been used to study the lattice 
dynamics of NissPd45, Ni5oPtso and NissCr^ alloys. In a more generalized context of 
inelastic neutron scattering in disordered alloys, an augmented space block recursion 
( ASBR) has also been proposed recently [9] . The ASBR calculates the full Green and 
Self-energy matrices instead of their diagonal entries alone. These are required for 
obtaining the response functions. 

In this paper, we shall make use of the ASBR technique to implement the 
theoretical formulation for the lattice thermal conductivity developed earlier by us. 
We shall study the dependence of lattice thermal conductivity on phonon frequency 
as well as on temperature in a series of disordered NiPd and NiPt alloys covering the 
full range of alloy compositions. 

As far as the temperature dependence of lattice conductivity is concerned, our 
results follow a general trend. We shall also discuss how the lattice conductivity 
behave as a function of concentration at several fixed temperatures. The notion of a 
minimum thermal conductivity will be justified. It will be shown that low temperature 
resonant modes considerably decrease the conductivity. An idea about the location 
of the mobility edge will be discussed from the phonon-frequency dependence of the 
mode averaged harmonic diffusivity. This dependence also gives a rough idea about 
what fraction of the states across the frequency spectrum are delocalized and therefore 
can carry current. The concentration dependence of the harmonic diffusivity will also 
be displayed. 

2. Theoretical formulation 

In the next two subsections, we shall apply the formalism introduced in our earlier 
paper (AM) to study NiPd and NiPt alloys across the alloy composition range. In 
(AM) we have discussed the theoretical formalism in great detail. We shall present 
only the main pointshere and leave the reader to refer to (AM) for greater detail. 

The derivation of a Kubo formula for thermal conductivity requires an additional 
statistical hypothesis, which states that a system in steady state has a space dependent 
local temperature T(r^) = [kbP^i)] -1 ■ The matrix element of the heat current in the 
basis of the eigenfunctions of the Hamiltonian is given by : 



here 7,7' label the various modes of vibration, z/k 7 ,fk 7 ' are their frequencies, 



s^ y (k) = - K + M < 7 ,(k) 



where, the phonon group velocity v 7 y (k) is given by 
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where 



e^(k), Cy (k) are the polarization vectors and D /ily (k) is the Fourier transform of mass 
scaled dynamical matrix. 

We shall consider the case where the temperature gradient is uniform within 
the system. The Kubo formula then relates the linear heat current response to the 
temperature gradient field 

/oo 
dt' ^{t-t') V v 5T{t), 
-oo 

^(r)=6(r) ijf" d\{S»(-ih\),S»(r)), 

9(r) is the Heaviside step function, and 

S(-ihX) = e XH S e- XH . 

( ) on the right-hand side of the above equation denotes thermal averaging over states 
in the absence of the temperature gradient. The above equation can be rewritten in 
the form of a Kubo-Greenwood expression 



7 7'#7 



d 3 k s-^ x ^ (nk 7 ') - (n k7 ) 



S^,(k)S 7 , 7 (k) 



^ [{/Sew s^(k)j {/ Sew s 77 (k)| 



^k v 9(n k7 ) 

8 ^3 2^ ^(/j^) 



5(1/) 



(1) 



(2) 



where (rik 7 ) = (eP hVk ~i — 1) 1 is the equilibrium Bose Einstein distribution function 
and T is the absolute temperature. 

The first expression is for inter-band transitions, while the second expression is for 
intra-band transitions. For an isotropic response, we can rewrite the first expression 
as 



Ki(l 



^' T ) = ^E j dv ' j 0E § 77'( k ' T ) S£ 7 (k,T)J(i/-i* y )5(i/ + I ,-i^) 



where 



S^(k,T) = 



(rik 7 ') - (n k7 ) 



/i(z/k 7 - ^k 7 ') 



S^(k). 
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KI ^ T ^ / dv ' I 5Tr[s^(k,r)3m{G(k,i,')}S^(k,r)3™{G(k,i,' + z,)} 

The operator G(v) is the phonon Green operator (Mv 2 l — <I>)~ 1 . The Trace 
is invariant in different representations. For crystalline systems, usually the Bloch 
basis {|k, 7)} is used. For disordered systems, prior to configuration averaging, it 
is more convenient to use the basis {|k, a)}, where k is the reciprocal vector and a 
represents the coordinate axes directions. We can transform from the mode basis 
to the coordinate basis by using the transformation matrices T 7Q (k) = e"(k). For 
example 



S^(k,T) = T- 7 T (k)S^(k,r)T^(k). 



If we define 



Tr[s G(k,zi) S G(k,z 2 )J . (3) 

then the above equation becomes, 

< v > T )=i^r £ / dv ' [« li,i ('/-y + + u) + K'">( V ' + y- + u) 

-K^{v' + ,V'+ + v)- ^{v'-y- + v)] (4) 



where 



/(!/+) = lim/(i/ + i<f), f{v-) = \\mf{v-i5). 

We have used the herglotz analytic property [16] of the Green operator 

G{v + it) = 3?e [GO)] - i sgn(J) 3?m [G(i/)] . 

For disordered materials, we shall be interested in obtaining the configuration 
averaged response functions. This will require the configuration averaging of quantities 
like k(zi,z 2 )- This averaging procedure we have discussed in detail in (AM) : 

< k(z 1 ,z 2 ) > = < «(i)(zi, z 2 ) > + < Ak(zi, z 2 ) laddcr > (5) 

The first term in the right hand side stands for the disorder induced corrections 
where the disorder scattering renormalizes the phonon propagators as well as the heat 
currents. In this term, corrections to the heat current is related to the self-energy of 
the propagators. The second term of Eq. (5) includes the vertex corrections due to 
the correlated propagation. 
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as 



D^(k) = 7T 4-S^,(k)S^(k) ^ ^ ^y). 

This is an intrinsic property of the 7-th normal mode and provides an unambiguous 
criterion for localization. 

The averaged thermal diffusivity (averaged over modes) is then given by 

d 3 k. 



7 7 

Assuming isotropy of the response, we can rewrite the numerator of above 
equation as 



where 



7 7 



1 



*-' (k) = ^ S -' (k) - 

We may again rewrite the above equation for D to t as 

1 f f r/ 3 k r ^ ~ 

ATO = y dv ' J ^3 Tr [^{G(ky)}S*(k)9m{G(ky)}S*(k)9to{G(k,i/)} 

The averaged thermal diffusivity can then be expressed as (for an isotropic 
response) 

dm -l£ d~m -\ . t ,3 k £ " D '°'"" «fl 

For disordered material, we shall be interested as before in obtaining the configuration 
averaged thermal diffusivity. The configuration averaged thermal diffusivity can be 
expressed (to a 1st order approximation) in the form 



J ^ Tr I 3m « G (k, v) »] 



where 



« DZ{v) » * ^\ dv ' \ S Tr [ 3m « G ( k ' ^ » 

< S"(k) 3?m {G(k,i/)} S"(k) 3ra{G(k,i/)} > 
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3. RESULTS AND DISCUSSION 

The details of numerical calculation for the two alloys of our interest in the present 
work are as follows : 

• We have carried out calculations on 501 ^-points. 

• A small imaginary part of the frequency 5 =0.001 has been used for evaluating 
the Green matrix and Self-energy matrix in the augmented space block recursion 



• The calculation of lattice conductivity has been done at 40 temperatures. 

• For the Brillouin zone integration, 145 k-points in the irreducible 1/48-th of the 
zone produced well converged results. 

3.1. NiPd alloy : Strong mass and weak force constant disorder. 

For a list of general properties of fee Ni and Pd, we refer the reader to Ref. [16]. 
This particular alloy has already been studied experimentally by Farrell and Greig [5] 
using conventional potentiometric techniques. But unfortunately their investigation 
was limited only to very dilute alloys in the temperature range 2-100 K. In an earlier 
communication (AM) we have already shown a comparison of our results for the 
temperature dependence of thermal conductivity with theirs for a dilute NiggPdoi 
alloy. 

Our initial focus in the earlier communication [1] was to calculate the effect of 
disorder scattering over an averaged medium. We had made use of the multiple 
scattering diagram technique and went beyond the framework of CPA which treats 
only diagonal disorder. Rather we had applied the technique in a more generalized 
context with the inclusion of diagonal as well as off-diagonal disorder arising out of 
the disorder in the dynamical matrix. In an attempt to compare the effect of averaged 
heat current and disorder induced renormalized heat current on the lattice thermal 
conductivity we have plotted Fig. 1. 

Fig 1 shows the results for disordered Ni^Pd^ alloy. The black curve represents 
the lattice conductivity including all kinds of disorder induced corrections : e.g. 
corrections to the heat current and the vertex corrections, while the red curve stands 
for the same quantity but using averaged heat currects and without vertex corrections. 
The green curve in this figure shows the scaled joint density of states. From the figure 
it is clear that the transition rate V is strongly dependent both on the initial and the 
final energies throughout the phonon frequency (v). That is 



[9]. 



K{y,T)^\r{y,T)\J{y) 



(8) 



where J{v) is the joint density of states given by 
J{v) = [ dv' [ P% Tr hm < G(k, v') »3m« G(k, v' + v) » 




(9) 



Figure 1 shows that, the effect of disorder corrections to the current terms on the 
overall shape of lattice conductivity is rather small in the Ni 50 Pd5o alloy. We should 
note that the effect of disorder corrections to current and the vertex corrections in 
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Figure 1. (Color Online) Configuration averaged lattice thermal conductivity vs 
phonon frequency v (THz) for NisoPdso disordered alloy. The red line and black 
line shows the conductivity using the average VCA current and effective current 
(consisting of average VCA current + disorder corrections + vertex correction) 
respectively. The green line in indicates the configuration averaged joint density 
of states. 



these alloys become negligible beyond phonon frequencies ~ 2.5 THz. For the higher 
frequency modes the effect of scattering phenomenon is well described by the mean 
filed approximation. It is in the low frequency region that configuration fluctuation 
effects beyond the mean-field becomes significant. 

There is a very important feature in Fig. 1 that still needs discussion : which is 
how to explain the origin of a dip in n,(y) at the lowest energy v = ? A similar kind 
of dip has also been reported by Feldman et.al. [17] while studying amorphous Si and 
Sii-^Gej. alloys. Their k(v) have a small Lorenzian shaped dip centred at i>=§. This 
reflects the missing intraband conductivity n 11 . This dip in n{v) stands at a small but 
finite value [v ~0]. The finiteness of the dip in n(v) is because of the fact that their 
calculation was based on a Kubo-Greenwood expression for the thermal conductivity 
with the delta functions in the expression broadened into a Lorenzian of small (but 
finite) width 77. However in our case it is evident from the Fig. 1 that this dip in k(v) 
stands at k(v) — > at v — 0. This is due to the simple reason that in our calculation 
the Lorenzian broadening has not been put in by hand, but it arises automatically 
from the disorder effect on the crystalline spectral function 3m[G(k,i/)] . Another 
reason for this difference in the position of dip in k{v) may be due to the fact that 
Feldman et.al. carried out their calculation at a fixed wavevector k, while we have 
summed over the entire Brillouin zone. 

The origin of this dip can also be explained by looking at the joint density of 
states (JDOS) represented by green line in Fig. 1. This quantity has a dip near v = 
reminiscent of the dip in the k(v) curves. This indicates that a smooth convolution 
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of two Green matrices G(k, v') and G(k, v' + v) (or two smooth densities of states 
obtained after summing k over the Brillouin zone), as appeared in the expression (5) 
of [1], is mainly responsible for such a sharp dip in the lattice conductivity at v = 0. 
As discussed by Feldmanei.a/, this dip at v — disappears as the system size N — > oo. 
They have also suggested an appropriate method to eliminate this dip in a sensible 
manner, which allows us to extrapolate the k(v) curve from a value at v > [y ~0) 
to a value at i/=0. This extrapolated value of n(v) at v = is nothing but the d.c. 
value of the lattice thermal conductivity kq. In an attempt to calculate (kq), we have 
obtained a value of 15.25 W/m/K for Ni 50 Pd 50 alloy at T=110 K. 




Figure 2. (Color Online) Thermal conductivity vs temperature T(K) for NiPd 
alloys and Amorphous Si. The top panel shows our results on the lattice 
conductivity for NiggPdoi alloy at three different frequency cut-off v. The middle 
panel shows the lattice conductivity for amorphous Si [17] at three different cut- 
off frequency, while the panel at the bottom shows the experimental data [5] for 
the total thermal conductivity (= lattice + electronic contribution) of the same 
NiggPdoi alloy. 



Direct comparison with the experimental data on these systems is difficult, 
because the experimental thermal conductivity also has a component arising out of 
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Figure 3. Residual or impurity contribution of the electronic part of the thermal 
conductivity 

the contribution from electrons. Figure 2 shows the temperature dependence of lattice 
conductivity. The top panel shows our theoretical result for the NiggPdoi alloy at 
three different frequencies. The bottom panel shows the experimental data [5] on the 
total ( electronic and lattice ) thermal conductivity of the same 99-01 NiPd alloy. Since 
the frequency is not mentioned in the experimental data, we assume that it must be for 
low frequencies. The best comparison then will be between the middle (black) curve 
on the top panel and that in the bottom one. The two agree qualitatively except at 
low temperatures where we expect the electronic contribution to dominate. In order 
to understand whether the deviation does arise from the electronic contribution, we 
have compared the top panel with the thermal conductivity of amorphous-Si [17], 
shown in the middle panel. In a-Si the electrons near the Fermi level are localized and 
hence cannot carry any current. The contribution to thermal conductivity arises from 
scattering due to configuration fluctuations in the amorphous material. Qualitatively 
we expect the results to be similar to configuration fluctuation scattering in random 
alloys. Almost the entire contribution should come from the phonons. The behaviour 
of the two panels are quite similar. The origin of the hump in the experimental 
lattice conductivity can also be understood if we assume Widemann-Franz law and 
write the residual part of the electronic contribution of the thermal resistivity as 
K r = LqT/ po(T). Here Lo is the Lorenz number and po(T) is the electrical resistivity. 
Assuming that the electrical resistivity behaves as po(T) = A + BT + CT 2 at low 
temperatures and with a suitable choice of the parameters, this contribution does 
show a hump followed by a decreasing behaviour flattening out at larger temperatures 
(see Fig 3). The sum of the contribution shown in the top panel of Fig 2 and that in 
Fig 3 would lead to the experimental behaviour shown in the bottom panel of Fig 2. 
This is a plausibility argument and needs to be confirmed by a detailed calculation of 
the electronic contribution to the thermal conductivity. 

We shall now present the dependence of lattice thermal conductivity and thermal 
diffusivity on parameters such as phonon frequency (v), concentration (x) etc. In Fig. 
4, we display the frequency dependence of lattice conductivity for Ni 50 Pd5o alloy at 
various temperatures evaluated using Eq. 5. The figure clearly shows the saturation 
of lattice conductivity as we proceed towards the higher temperatures. The d.c. value 
of the conductivity (kq), which is just extrapolation of k(u) curve from a value at 
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Figure 4. (Color Online) The configuration averaged lattice thermal conductivity 
vs phonon frequency v (THz) at different temperatures T for NisrjPdso alloy. 
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Figure 5. (Color Online) The averaged lattice thermal conductivity vs 
temperature T(K) at various cut-off frequencies v cu toff f° r NisoPdso alloy. 



v > to a value at v = 0, increases as we increase the temperature. 

The temperature dependence of lattice thermal conductivity for the NisoPdso 
alloy at various phonon frequencies v are shown in Fig. 5. It is qualitatively similar 
to NisoPtso alloy. The conductivity initially increases with temperature reflecting a 
quadratic T-dependence (in the low T-regime) and ultimately reaches a T-independent 
saturated value (at higher temperatures). 



Figure 6 shows the lattice conductivity as a function of frequency at T=100 K 
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Figure 7. (Color Online) Lattice thermal conductivity vs Pd-concentration for 
various temperature T at phonon frequency v = l.ObTHz. 

for various alloy compositions. Comparing the results of Fig. 6 with those of Fig. 14 
for NiPt, it is clear at a glance that the overall shape of frequency dependence of k for 
various alloys of Ni 1 _ 2 .Pd 2; looks similar, however for x=0.8 and a;=0.9 extra structure 
appears in the frequency dependence. Similar behaviour has also been observed for 
ir=0.9 in Nii-^Ptx alloy in Fig. 14. We believe that this behaviour for Nii-^Pdx 
alloy at x=0.8 and 0.9 may be due to the strong disorder in masses, the effect of 
which becomes important in the two dilute limit alloys. Such an anomalous behaviour 
is also reflected in the concentration dependence of lattice conductivity. This is shown 
in Fig. 7, which plots the lattice conductivity vs Pd-concentration (x) at a fixed 
frequency ^=1.05 THz for various temperatures. 




Figure 8. The configuration averaged thermal diffusivities D(v) for Nii-^Pd^ 
alloys, (a) x=0.1; (b) x=0.3; (c) x=0.5; (d) x=0.7; (e) x=0.9.The broad line on 
the frequency axis shows the extent of the vibrational spectrum. 



The thermal diffusivities D(^) are important because the effect of disorder is often 
manifested in them more directly than in the conductivities. Not only that thermal 
diffusivity also gives an approximate idea about the location of mobility edge as well as 
the fraction of delocalized states. In Fig. 8 we display the thermal diffusivity D{v) vs 
frequency for various compositions. The extent of the phonon frequency spectrum is 
shown by the broad lines. The density of states is non-zero across this spectrum range. 
The first thing to note is that, the region of large diffusivity in the five sets of alloys at 
the higher frequency side is not the same. In other words the weakly defined hump in 
the lower as well as higher frequency side are located at different positions for different 
alloys. The low frequency maximum in diffusivity is a minimum around the 50-50 
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Figure 9. The position of the mobility edge (bottom) and the percentage of 
mobile phonon states (top) as a function of the alloy composition for NiPd alloy. 



composition where the disorder scattering is the maximum. Above 2.8 THz, there is 
a smooth decrease of diffusivity approximately linear in frequency D(^) oc (y c — v) a 
, with the critical exponent a ~ 1 and a critical frequency v c where D{v) vanishes to 
within a very small level of noise. The allowed phonon states beyond this frequency 
must be due to localized phonon modes. The critical frequency v c locates the mobility 
edge above which the diffusivity is strictly zero in the infinite size limit. Once the 
mobility edge is located, the fraction of de-localized states may be estimated directly. 
It is clear from the figure that location of mobility edge varies with composition. 
Consequently the percentage of de-localized states available for thermal conduction 
in the system also varies with composition. An inspection of Fig. 8 determines the 
location of the mobility edges (y c ) for the five different compositions. The following 
figure 9 which shows the position of the mobility edge and the percentage of mobile 
phonon states in the spectrum as a function of the composition is quite illustrative. 

The maximum percentage of localized states occur at 50-50 composition where 
we expect disorder scattering to be a maximum. The mobility edge moves to higher 
frequencies as the concentration of Ni increases, but so does the band width of the 
phonon spectrum. The minimum percentage of mobile phonon states available for 
thermal conduction occurs, as expected, at around the 50-50 composition. A similar 
behaviour has also been discussed by Feldman et.al [17] while studying the effects of 
mass disorder on various Sii^Ge^ alloys. 
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Figure 10. (Color Online) The averaged thermal diffusivity D(u) vs Pd- 
concentration at various cut-off frequencies v cuto jj for Nii—^Pdz alloy. 



Figure 10 shows the concentration dependence of thermal diffusivity for Nii^^Pd^ 
alloy at various frequencies (v). As expected from our earlier discussion, the 
minumum diffusivity occurs around the 50-50 composition, where disorder scattering 
is maximum. The curves have asymmetry around x — 0.5 which decreases with 
increasing frequency. This asymmetry reflects a similar asymmetry in the thermal 
conductivity as a function of composition. 



3. 2. NiPt alloy : Strong mass and force constant disorder. 

We refer the reader to a previous article [16] by us for some of the basic properties of 
fee Ni and Pt which is relevant for our present calculation. It has been our experience 
[16, 9] that the effect of disorder in NiPt alloy is more dramatic than NiPd. For 
instance the appearance of sharp discontinuities observed in the dispersion where we 
have resonance states and consequent increase in the line-width [16]. 

Figure 11 shows the results for disordered NisoPtsa alloy. As before, the 
black curve represents the lattice conductivity including all kinds of disorder induced 
corrections : e.g. corrections to the heat current and the vertex corrections, while the 
red curve stands for the same quantity but using averaged heat currents and without 
vertex corrections. The green curve shows the scaled joint density of states. From the 
figure it is clear that as in the case of NiPd, the transition rate V is strongly dependent 
both on the initial and the final energies throughout the phonon frequency (y) . Figure 
11 also clarifies that although the effect of disorder corrections to the current terms is 
small, but this effect is comparatively more pronounced in Ni^Pt^ than in Ni^Pd^Q. 
The effect of disorder corrections to current and the vertex corrections in the present 
case become negligible beyond phonon frequencies ~ 2.4 THz. In the low frequency 
region, configuration fluctuation effects beyond the mean-field is more pronounced in 
NiPt than in NiPd. This may be because of the two simple physical reasons : First, the 
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Figure 11. (Color Oncurve) Configuration averaged lattice thermal conductivity 
vs phonon frequency v (THz) for NisoPtso disordered alloy. The red curve 
and black curve shows the conductivity using the average VCA current and 
effective current (consisting of average VCA current + disorder corrections+ 
vertex correction) respectively. The green curve indicates the configuration 
averaged joint density of states. 



NiPt is an alloy where both mass (mp ( / rxijvi ~ 3) as well as force constant (Pt-force 
constants are on an average 55% larger than those of Ni) disorder dominates, while 
in NiPd alloy the mass disorder (rnpa/ rnjv, ~ 1.812) is weaker than NiPt and the 
force constants are almost the same for the two constituents. Second, from a purely 
phenomenological point of view, there is a larger size mismatch between Ni and Pt in 
NiPt alloy as compared to Ni and Pd in NiPd alloy. 

In Fig. 12 we display the frequency dependence of lattice thermal conductivity 
for Ni5oPt5o alloy at various temperatures, evaluated using the formulae given in Eq. 
(3) and (5). We note from Fig. 12 that as we increase the temperature, the difference 
in lattice conductivity at a particular frequency is quite large in the lower temperature 
regime ( 20 < T < 80K) , however the difference starts to saturate as we go to higher 
temperatures. As discussed in the previous section, the extrapolated value of k{v) 
curve from a value at v > Q ( next to v = ) to a value at v = is called the 
d.c. value of lattice thermal conductivity (kq). Figure 12 clearly shows that this d.c. 
value of the conductivity increases as we go on to increase the temperature. This is 
expected because in the d.c. limit the mechanism of heat conduction is transfer of 
energy between delocalized modes of equal energy by the heat current operator. 

The behaviour of k(u) vs v curve in both NiPd and NiPt is qualitatively similar, 
but there are some quantitatively different features. The maximum in NiPt is located 
at a lower frequency as compared with NiPd. Also the overall magntitude of lattice 
conductivity in NiPd is higher than NiPt. 

The temperature dependence of lattice conductivity for NisoPtso alloy at various 
cut-off frequencies are shown in Fig. 13. The conductivity increases initially ( in the 
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Frequency V (THz) 



Figure 12. (Color Online) The 
conductivity vs phonon frequency v 
NisoPtso alloy. 



configuration averaged lattice thermal 
(THz) at different temperatures T for 



low T-regime ) as an approximate quadratic function of temperature and ultimately 
increases smoothly to a T-independent saturated value. As far as such dependence of 
k(T) in the high temperature regime is concerned, the heat in this conduction channel 
is carried by non-propagating modes which are strongly influenced by the disorder 
but mostly not localized and therefore able to conduct by intrinsic harmonic diffusion. 
This is a smooth dependence which closely resembles the specific heat and saturates 
like the specific heat at high temperaturs. Following Slack [4] we could call this piece 
the "minimum thermal conductivity . A number of authors have discussed that low 
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Figure 13. (Color Online) The averaged lattice thermal conductivity vs 
temperature T(K) at various cut-off frequency v cuto jj for NisoPtso alloy. 
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Figure 14. (Color Online) Frequency dependence of lattice thermal conductivity 
for various alloys Nii—^Pta; at T=100 K. 

temperature dependence of k(T) shows a mild plateau like region. In this regime 
the heat is mainly carried by the propagating long wavelength acoustic modes. The 
complex inelastic scattering processes then kill off the low frequency contribution at 
higher temperatures leaving a peak which becomes the plateau. However there are 
situations, where the propagating modes become reasonably well damped ( as in our 
calculation ) and are no longer able to carry much heat. In such cases the contribution 
of delocalized and poorly conducting vibrations takes over, giving a net result in good 
accord with the Kittel's old idea. Under these circumstances the plateau like region in 
the low temperature regime almost disappears. The damping of propagating modes 
is also amplified as we make the alloy more and more concentrated. This can easily 
be verified by looking at the results of reference [5] , which shows that as we increase 
the concentration the plateau like region goes down and gets smoother. 

In Fig. 14, we display the frequency dependence of lattice conductivity for NiPt 
alloys at various compositions, but at a fixed temperature T=100 K. As before each 
of the curves have a dip at the lowest frequency. It is important to notice that as 
we increase the Pt-concentration x, structure appears in the behaviour of k(v). This 
indicates that the structure arises due to the contribution of Pt-atom (with large mass) 
in the alloy. 

The concentration dependence of lattice conductivity at a fixed phonon frequency 
v = l.ObTHz are plotted in Fig. 15. The various curves in this figure stand for various 
values of the temperature T starting from a lower value of 50 K to a higher value of 230 
K. It is clear from the figure that the concentration dependence is almost symmetric 
about x=0.5. It has been discussed by Flicker and Leath [14] within the framework 
of coherent potential approximation that this asymmetry is a function of the size of 
the sample chosen i.e. a large N leads to less asymmetry. They have verified this 
statement by performing two calculations one for N=100 and other for N=2000. The 
concentration dependence in the later case is more symmetric as compared to the 
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Figure 15. (Color Online) Lattice thermal conductivity vs Pt-concentration for 
various temperature T at phonon frequency v = l.ObTHz. 



former ones. In our case, the results shown in Fig. 15 are the optimal symmetric 
structure for the concentration dependence of k. This is because our calculations are 
performed in the reciprocal space representation which involves the entire lattice. 

An interesting challenge remaining in this problem is to calculate the effect of 
adding anharmonicity to the model. The reason why one should be interested in 
calculating the anharmonicity effect is because in real systems at high temperature 
phonon-phonon Umklapp scattering becomes the dominant scattering mechanism. 
This Umklapp scattering actually arises due to the presence of anharmonic terms in 
the Hamiltonian. The effect of this anharmonicity is to flatten the lattice conductivity 
vs concentration curve. 

Fig. 16 shows the frequency dependence of diffusivity D{v) for various alloy 
compositions. The thick lines on the frequency axes shows the extent of the frequency 
spectrum. These have been obtained from the density of states calculations presented 
in our first paper on the subject [16]. It is clear from the Fig. 16 that there are 
basically two regions of large thermal diffusivity : one near the lower frequency region 
(~ Q.bTHz) and the other around a somewhat higher frequency region (~ \.2bTHz). 
But as we go on increasing the Pt-concentration, the former region of large diffusivity 
starts decreasing gradually and becomes almost flat for the maximum Pt-concentration 
of 90%. The latter region of large diffusivity sits on a portion of the frequency spectrum 
above the transverse acoustic vibrations. Here the modes have large velocities and 
are probably very effective carriers of heat. The approximate linear decrease in 
diffusivity starts at ~ 3 THz. But the location of mobility edge in this case varies 
with composition in a slightly different way as compared to the case of NiPd alloy. 
Fig. 17 (bottom) shows the position of the mobility edge v c as a function of the 
alloy composition. As the concentration of heavy Pt increases, the band width of the 
frequency spectrum (which is proportional to the square root of the mass) shrinks and 
the position of the mobility edge within the band also shrinks. Fig. 17 (top) shows the 
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Figure 16. The configuration averaged thermal diffusivities D(v) for Nii—^Pta; 
alloys, (a) x=0.1; (b) x=0.3; (c) x=0.5; (d) x=0.7; (c) x=0.9 .The broad line on 
the frequency axis shows the extent of the vibrational spectrum. 

fraction of the frequency band which is extended. When the disorder is the strongest, 
i.e. at 50-50 composition, this fraction is a minimum. 

One thing is very clear from the above discussion that, in an alloy where mass 
disorder dominates and the forceconstant disorder is weak ( as in the case of NiPd alloy 
), the complex disorder scattering processes try to localize more vibrational modes as 
compared to those in an alloy where both mass as well as the force constant disorder 
dominates ( as in the case of NiPt alloy ) . The result can be interpreted in a slightly 
different way as : "the role of force constant disorder in binary alloys is to make the 
vibrational eigenstates more delocalized " , i.e. the more dominant the force constant 
disorder is, the more delocalized the vibrational modes will be. 
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Figure 17. The position of the mobility edge (bottom) and the percentage of 
mobile phonon states (top) as a function of the alloy composition for NiPt alloy. 



In Fig. 18, we present the concentration dependence of thermal diffusivity for 
various phonon frequencies v. As in NiPd, the diffusivity are asymmetric about 
x=0.5. However this asymmetry reduces as we increase the phonon frequency. This 
asymmetry is reflected also in the behaviour of thermal conductivity with alloy 
composition. 
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Figure 18. (Color Online) The averaged thermal diffusivity D(v) vs 
concentration at various cut-off frequencies u cuto ff for Nii—^Ptx alloy. 
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4. Conclusions 

We have performed a detailed numerical study of the theoretical formulation (AM) 
developed earlier by us for the lattice thermal conductivity of disordered binary alloys. 
We have demonstrated through our numerical results that how this multiple scattering 
based formalism captures the effect of off-diagonal and environmental disorder present 
in the problem. The use of augmented space method to keep track of the configuration 
of the system and the block recursion method have made the implementation simple 
yet powerful. This is reflected in the satisfaction of essential herglotz analytic property 
of the diagonal green function in our earlier calculation [16]. A significant contribution 
of this article beyond the earlier theoretical approaches is the inclusion of force 
constant fluctuations properly in the theory. We have applied the formalism (AM) to 
two real disordered alloys ; namely NiPd and NiPt. We have shown that the effect 
of disorder corrections to the current and the vertex correction on the overall shape 
of lattice thermal conductivity for both the alloys are very small. Comparatively the 
effect is found to be more pronounced in NiPt alloy, which is due to the presence of 
strong disorder both in masses and force constants in this alloy. The prominence of 
the force constant disorder in NiPt alloy has also been demonstrated in the frequency 
dependence of lattice conductivity for various compositions of the Nii-^Pt^ alloy. 
The saturation of lattice conductivity at higher temperatures has been shown for 
both the alloys. The concentration dependence of k in Nii-^Pd^ alloy has been 
shown to be more asymmetric about x=0.5 than in Nii-^Pt^ alloy. The numerical 
results on the harmonic diffusivity provide an interesting idea about the localization 
and derealization of the vibrational eigenstates. It says that in disordered binary 
alloys " the more stronger the force constant disorder is, the more delocalized the 
vibrational modes will be" . That is why NiPt alloy has larger fraction of delocalized 
states as compared to that in NiPd alloy. For both the alloys, however we had no prior 
information about the species dependence of the force constants but rather choose a 
set of force constants intuitively as we have done earlier [9,16]. A better understanding 
of the role of disorder in the transport properties of random alloys could be achieved 
with prior information about the force constants. These could be obtained from more 
microscopic theories, e.g., the first principles calculation on a set of ordered alloys. 
Our future endeavor would be to rectify this and attempt to obtain the dynamical 
matrix itself from such microscopic theories. 
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